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1. Introduction 


Air pollutants are one of the causative agents in the presently seen forest decline syndrome, but 
the mechanisms of action are incompletely known. The classical acid rain hypothesis states that 
deposition of air pollutants contributes to soil acidification, loss of plant nutrients and mobilization 
of aluminium (ULRICH, 1983). In various areas, acidification of soils induced by atmospheric 
deposition indeed surpasses the internal rate of proton production (VAN BREEMEN etal., 1984). This 
soil-centred theory is often presented as an alternative to the direct damage caused by the action of 
reactive air pollution components on the leaves (ARNDT et al., 1982). Since diagnostic fertilization 
may reduce the symptoms of forest decline (HUTTL, 1985), it is clear that nutrient shortage, either 
induced by leaf damage, or by soil impoverishment, is of particular importance. To discriminate 
between the various theories of forest decline, one might argue that, if forest decline is at least partly 
due to processes occurring in the soil, the deterioration will be accompanied by changes in the soil 
microfauna composition, because of the intimate contact of soil animals with the soil. In this paper 
we investigate the relation between forest vitality, soil chemistry and the density of microarth- 
ropods. 

A number of authors have described changes in the soil fauna community as a consequence of 
artificially applied acid rain (BAATH et al., 1980; HAGvAR & AMUNDSEN, 1981; HAGvar, 1984; 
HEUNGENS & VAN DAELE, 1984). These studies demonstrate that rearrangement of the community, 
induced by soil acidification, may be very complicated. Among Collembola, the species Isotoma 
notabilis appeared to be sensitive to acidification, and could possibly serve as an indicator species. 
However, some of the results of these studies are in apparent conflict with each other (cf. HAGVAR, 
1984; HEUNGENS & VAN DaELE, 1984). This could be due to circumstances specific for the studied 
sites, or to the rate of acidification under artificial acid rain. 

In this paper we try to resolve the problem by means of a comparative approach, thus avoiding 
possible effects specific to one site, and effects from artificial acidification. We selected various 
forest stands of differing vitality, within an area of restricted size, which are about equally exposed 
to air pollution. To maximize the applicability of the approach, we focussed attention on easily 
recognizable arthropods, and expressed their population density as a percentage of the total density 
of the group to which they belonged; as argued by CLAUSEN (1986), relative population densities 
may be better indicators of pollution than absolute densities. On the basis of the results, we argue 
that monitoring of sensitive soil microarthropods (indicator species), in connection with chemical 
analysis, may contribute to judging the probable vitality development of a pine stand in a relatively 
early stage. 
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2. Material and methods 


2.1. Site selection 


In May 1985, 12 stands of Pinus sylvestris L. were selected, out of 22 inspected, within the Veluwe area in the 
Netherlands. The stands differed in vitality but were otherwise more or less equal with respect to age (> 20 years), soil 
type (podzol developed on sand), and ground vegetation (mainly consisting of Avellana flexuosa TRIN., Vaccinium 


Fig. 1. Map of the Veluwe area, showing the locations of sampling sites. The inset shows the Netherlands. 


myrtillus L., Dicranum scoparium HEDW. and Hypnum cupressiforme L.). Vitality of a forest was judged on the basis 
of the monitoring program of the Netherlands State Forest Management (Anonymous, 1984), as well as on own 
observations. Fig. 1 shows a map of the locations. . 

In each stand, a representative plot (20m Xx 20m) was delineated, and 15 random core samples were taken in the 
plot, by means of a split-corer (Ø 10cm). Sampling was confined to the litter-fermentation layer of the soil profile 
(Ago-Ao horizons combined) ; the depth of this layer was measured for each core, after opening the corer. For each plot, 
3 cores were used for soil chemistry aspects and 12 were extracted for microarthropods (see below). Vitality of the 
plot was measured as the number of year-classes of green needles in the crown of each of 10 randomly selected trees, 
using the scale: 1/2, 3/4, ... 2%. Each tree was judged by 3 people from different points of view. 
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2.2. Soil chemistry 


Each core was homogenized and pH was measured in a water extract of the fresh material (20g sample + 200ml 
deionized water, shaken for 2h). From each core, 2 X 100mg of oven-dried material was digested completely in a 
HNO,/HCIO, mixture (7: 1). Concentrations of K, Ca, Mg, Fe and Mn in the diluted digests were determined by flame 
atomic absorption spectrometry (Perkin Elmer 4000). Aluminium was determined in the digests by means of a 
colorimetric method (methylthymolblue, following reduction of Fe by ascorbic acid). Total amounts of carbon and 
nitrogen were determined by elemental analysis on | mg of oven-dried material (Carlo Erba 1106). On each of the 
3 replicate cores from a plot, determinations were done in duplicate, which, in the data analysis, were combined to a 
single observation for a core. 


2.3. Soil fauna 


Arthropods were extracted from the cores by means of Tullgren-type equipment (VAN STRAALEN & RUNINKS, 
1982). The extracted fauna was identified to groups recognizable under a dissecting microscope (magnification 
20—40 x). Total numbers of the orders Collembola and Oribatida were recorded for each core, as well as the numbers 
of the families Pelopidae, Damaeidae, Camisiidae, Phthiracaridae, Onychiuridae, Poduridae, Isotomidae and 
Sminthuridae. Platynothrus peltifer Kocu (Oribatida) and [sotoma notabilis SCHAFFER (Collembola) were recorded at 
the species level. Identification followed SELLNICK (1960) and Gistn (1960). Counts of microarthropods were 
expressed as numbers per core, and as a percentage of the total count of oribatids or Collembola in the same core. 


2.4. Data processing 


The following rationale was applied in searching for patterns in the data: 


1. Each variable was subjected to one-way analysis of variance, to test for differences between sites, regarding 
differences between cores within a site as sampling error. If necessary to achieve homogeneous variance, 
transformations were applied to the data. The logarithmic transformation was used for counts, and the arcsine 
transformation for percentages (SOKAL & ROHLF, 1981). 

2. For variables showing significant differences between sites, a pooled value was calculated for each site. If a 
variable had been transformed, the pooled value was calculated from the untransformed data, and was then 
transformed. These averages (12 for each variable) were tested for correlations with each other (Pearson’s test, 
twotailed). = 

3. From the variables showing one or more significant correlations with other variables, a cluster of variables was 
selected, which were correlated with each other, as well as with vitality. On this cluster the main conclusions were 
based. 


3. Results 


Vitality of the various tree stands, as judged by the number of needle year-classes in pine trees, 
is given in table 1. Taking the introduced genotypes and local factors into account, the presence of 
two year-classes is considered normal for Pinus sylvestris stands in the area (Anonymous, 1984). 
Our survey was made in early spring, before budding of that year’s needles. In late summer, the 
number of year-classes is about 1/2 higher (unpublished observations). Given these conditions, it 


Table 1. List of the investigated sites, their codes, and the vitality of their tree stands, as judged by the mean number 
of needle year-classes in May 1985, with standard deviations based on 10 trees per site 


site code local name of forest needle 
year-classes + sd 


A Landgoed Welna 1.35 + 0.41 a 
B Boswachterij Nunspeet 0.85 + 0.34 b 
C Willemsbos 0.93 + 0.24 b,c 
D Leuvenhorst 1.10 + 0.21 a,c 
E Leuvenumse bos 1.00 + 0.00 b,c 
F Sprielder bos 1.25 £0.35 a,c 
G Boswachterij Kootwijk 0.90 + 0.21 b 
H Loobos 0.85 + 0.23 b 
I Westerrode 0.88 + 0.18 b 
J Boswachterij Hoenderloo 1.03 + 0.30 b,c 
K Elspeter bos 0.88 + 0.21 b 
L Stakenberg 0.85 + 0.21 b 


Note: Means sharing a common letter (a, b, c) do not differ significantly at the 5% level. 
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appears from table 1 that one half of the investigated pine stands have low vitality (sites L, H, B, 
K, I and G). In fact only the stands at sites (A, F and D) can be considered as clearly vital, while 
those at sites E, J and C take an intermediate position. 

Table 2 provides a survey of soil variables which did not differ significantly between the sites, 
while table 3 lists the variables that did show significant differences. Comparing these tables, it 
appears that the sites form a rather homogeneous set with respect to the elements Ca, Mg and Fe, 
but are different in the elements K and Mn. The pH does not differ: all soils are equally acid. The 


Table 2. List of soil variables not differing significantly between sites 


variable mean + sd units 
depth litter layer 622 13 cm 

pH (H20) 37 + 0.1 ra 

Ca 30.7+ 9.1 umol/g 
Mg 12.2 = 22 umol/g 
Fe 52.8 + 11.6 uumol/g 
Al 16.5+ 2.2 umol/g 


Note: Means and standard deviations were based on observations from all sites combined. All variables refer to the 
A-horizon of the forest floor. 


nitrogen content of the forest floor was highest for the sites A, I, F and H. Since this includes both 
vital (A, F) and less vital (I, H) sites, there is no clear relation of forest floor nitrogen to vitality. A 
similar argument holds for K. 


Table 3. List of soil variables showing significant differences between sites 


variable range units 

K 4.6— 9.9 umol/g 

Mn 0.73= 2.6 umol/g 
— N 0.83— 1.44 % (m/m) 

C/N 22—30 = 


Note: For each variable a range is given (site with lowest mean — site with highest mean). All variables refer to the 
A-horizon of the forest floor. 


Among the variables given in table 3, the Mn concentration and the C/N ratio were 
significantly correlated with vitality from table 1. This is pictured in figs. 2 and 3. The most vital 
sites (A and F) have the highest Mn content and the lowest C/N-ratio, while the other sites are 
more like each other. 

Among the soil fauna, almost all groups recorded showed significant differences between the 
sites. The densities of total Oribatida and Collembola vary from 8,000 to 80,000 per m? 
(Oribatida) and from 4,000 to 30,000 per m? (Collembola). There is a clear correlation between 
Oribatida and Collembola (fig. 4), but the density ratio (on the average 3.6:1) was not 
significantly correlated with vitality. It seems that the sites provide habitats which support 
different densities of microarthropods, for reasons not always related to the vitality of the stand. 
For some of the groups (notably Collembola) the depth of the litter layer could explain some of the 
variation. To avoid aspecific effects of site differences, we consequently expressed the density of 
each arthropod group as a percentage of the density of its order (Oribatida or Collembola). 

From the various density ratios tested, two were significantly correlated with vitality: the 
number of Platynothrus peltifer as a percentage of the Oribatida, and the number of Sminthuridae 
as a percentage of the Collembola. In the two vital sites A and F, P. peltifer accounts for 10—20 % 
of the density of oribatids, while in the other sites this percentage was only about 2% (fig. 5). 
Other oribatid groups were more abundant at these sites, but not in a consistent pattern related to 
vitality. Among the Collembola, the proportion taken by the family Sminthuridae varied from 9 to 
65% and was correlated with vitality (fig. 6). The other main collembolan group was Isotomidae 
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Fig. 2. Correlation between vitality and the Mn concentration of the forest floor. Each point refers to the mean of a 
site. For site codes see Table 1. Correlation coefficient = + 0.74. 
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Fig. 3. Correlation between vitality and the C/N ratio of the forest floor. Each point refers to the mean of a site. For 
site codes see Table 1. Correlation coefficient = — 0.62. 


and this family showed the opposite trend. The species /sotoma notabilis was present in low density 
at most of the sites and could therefore not serve as an indicator. 

Comparing the results on microarthropods and soil chemistry, it can be concluded that, although 
the correlations are weak in some cases, there is a consistent pattern in the data: the soils of the most 
vital stands are highest in Mn, lowest in C/N, and have the highest relative population densities of P. 
peltifer and Sminthuridae. 


4. Discussion 


The analysis of the relations between soil microarthropods and soil fertility has a long tradition, 
but the application of this study to the present-day problems of acid precipitation and forest decline 
adds a new dimension to it. Our results have shown that, within the study area, both soil fauna and 
soil chemistry aspects differ between pine stands of differing vitality in a consistent manner. 
Therefore, we conclude that changes in the soil are at least part of the forest decline problem in the 
Veluwe area in the Netherlands. 

Due to the inventory nature of the study, precise relations between forest vitality, soil quality 
and soil fauna, are still incompletely known. Among the soil chemistry aspects, the carbon/nitrogen 
ratio and the Mn concentration were correlated with vitality and these will be discussed first. 
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Fig. 4. Correlation between the mean population densities of Collembola and Oribatida. For site codes see Table 1. 
Correlation coefficient = + 0.78 (both axes logarithmically transformed). 
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Fig. 5. Correlation between vitality and the population density of Platynothrus peltifer, expressed as a percentage of 
the total number of Oribatida. For site codes see Table 1. Correlation coefficient = + 0.75; y-axis transformed 
according to y = arcsin (V(y/100)). 
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Fig. 6. Correlation between vitality and the population density of Sminthuridae, expressed as a percentage of the total 
number of Collembola. For site codes see Table 1. Correlation coefficient = + 0.70 (y-axis transformed as in 
Fig. 5). 
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The higher average C/N-ratio as found in the forest floor of less vital stands, may relate to a 
disturbance of litter decomposition at these sites. It is well known that nitrogen is conserved during 
the initial stages of decomposition, and that the C/N-ratio tends to decrease (SwiFT et al., 1979). A 
high C/N ratio is therefore indicative of an over-representation of undecomposed material in the 
litter. Acid precipitation may inhibit microbial activity in soil (BEWLEY & StoTzxKy, 1983), and 
this may affect both immobilization and mineralization of nitrogen, depending on the stage of 
decomposition of the litter (BosaTTaA, 1982; KLEIN et al., 1984). Since soil microarthropods take 
part in the nitrogen dynamics of pine needle litter by grazing microbial biomass (VERHOEF & DE 
GoEDE, 1985), differences in microbial activity will be parallelled by changes in the microarth- 
ropod community. Another disturbance of the nitrogen dynamics may be the ammonia emission 
from pig sties, chicken farms and field-spread dung, which is rather intense in the area. Pine trees 
exposed to deposition from these sources have an increased nitrogen content in their needles; also, 
litter from these places shows a disturbed pattern of decomposition and nitrogen dynamics 
(VERHOEF & DOREL, 1988). Effects from ammonium deposition have recently been considered as 
an additional explanation of forest decline (NIHLGARD, 1985). In our study, these emissions have 
probably affected all sites in an equal manner, since the non-diffuse ammonium deposition from 
agricultural sources does not exceed 1,200m (ROELOFs et al., 1985), and none of our sites was 
contiguous to a farmyard or field. 

As to the manganese concentrations, the differences between sites probably reflect small 
differences in soil fertility. Being a relatively mobile element, Mn is one of the first nutrients 
(preceded by K and Mg) to leach under acidification (HOVLAND et al., 1980; BROWN, 1985). Under 
acid conditions, there is often a net loss of Mn from a forest ecosystem, because the output by 
leaching is not compensated by an equal input from precipitation (ZOTTL, 1985). It is therefore 
conceivable that vitality differences between sites are partly due to nutrient shortage in the soil. 
Litter Mn concentrations may also be influenced by the fact that Mn is among the elements 
accumulating in pine needles by aging (ERNST, 1985), so that a premature loss of needles may lead 
to a decreased input of Mn to the forest floor. 

Soil acidification, in the sense of depletion of acid neutralizing capacity, may proceed for a long 
time without change in soil pH (VAN BREEMEN et al., 1983); therefore, it is not surprising that soil 
pH in the present study did not differ between sites, although their acid neutralizing capacity may 
differ due to small differences in soil type. How soil acidity varies on a scale comparable to the 
microhabitat of soil invertebrates remains unknown. 

Acidification of forest soils by simulated acid precipitation may lead to both negative and 
positive effects on various soil invertebrates (BAATH et al., 1980; HAGvAR & AMUNDSEN, 1981; 
ABRAHAMSEN, 1983; HUHTA, 1984; HÅGVAR, 1984; HEUNGENS & VAN DAELE, 1984). Also, 
invertebrates may suffer directly from gaseous air pollution, although there appear to be large 
differences between species (LEBRUN ef al., 1978; LEETHAM et al., 1982; ANDRE et al., 1982; 
WHITE, 1983; WEIGMANN & Kratz, 1987). According to HAGVaR (1984), anthropogenically 
induced acidification of forest soils leads to community changes in favour of species with an already 
high dominance in naturally acid soils. The collembolan Jsotoma notabilis seemed to be relatively 
acid-avoiding, as judged from its distribution in Norwegian forest soils (HAGVAR & ABRAHAMSEN, 
1984), and this was consistent with its decrease under experimental acidification (HAGVAR, 1984). 
In our study area, hardly any Isotoma notabilis was present, although this species is abundant at 
various other places in the Netherlands ; this may relate to the rather extreme situation in the Veluwe 
area with an average litter pH (H20) of 3.7. It is unknown whether sminthurids, which were found to 
correlate with vitality in this study, are more acid-sensitive than other Collembola. This family is a 
group of rather active, superficially living springtails, which, because of their drought-resistance, 
would seem to be less dependent on taking up free water than other Collembola. Therefore, a direct 
effect of soil pH does not seem probable for this group. Rather, results obtained in 1986 (VAN 
STRAALEN ef al., 1987) suggest that, with decreasing vitality there is an increase of total 
Collembola, without an increase of Sminthuridae. It would seem that other collembolan groups 
(Isotomidae, Entomobryidae) profit by acidification and thereby cause the proportion of sminth- 
urids to decrease. 

For the oribatid mite Platynothrus peltifer, experimental research (VAN STRAALEN et al., 1987), 
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has shown that this species avoids acid substrates. Also, P. peltifer appears to have a markedly high 
body content of Fe and Mn, compared to Collembola: ithas 10 times more Fe and 50 times more Mn 
than the collembolan Orchesella cincta [VAN STRAALEN etal., 1987]. It would seem that, because of 
its acid sensitivity, and its high need for micronutrients such as Mn, P. peltifer is on the edge of its 
ecological amplitude in the Veluwe area, and is decreasing in abundance at those sites where the forest 
soil is critically depleted of nutrients and where preferred (relatively alkaline) microsites are 
becoming scarce. It is remarkable that the sites which take an intermediate position with respect to 
vitality (J, E and C) resemble the less vital sites more than the vital sites, when considered with respect 
to P. peltifer (fig. 5). Therefore, we conclude that this mite could be a valuable indicator species for 
monitoring changes in soil quality related to forest decline in an early stage. 
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Synopsis: Original scientific paper 
STRAALEN, N. M. van, M. H. S. Kraak & C. A. J. DENNEMAN, 1988. Soil microarthropods as indicators of soil 

acidification and forest decline in the Veluwe area, the Netherlands. Pedobiologia 32, 47—55. 

A comparative study has been made of the abundance of various soil microarthropod groups in 12 coniferous 
forests within the Veluwe area, the Netherlands. Also, a number of nutrient elements in the A-horizon of the forest soil 
have been measured. Relationships between forest vitality (estimated by number of needle year-classes), abundance of 
microarthropods, and soil nutrient concentrations were studied by correlation analysis. The results show that the 
oribatid mite Platynothrus peltifer and the collembolan family Sminthuridae are most abundant (relative to total 
Oribatida and total Collembola) at the most vital sites and decrease as forest vitality decreases. This is correlated with a 
decrease of the Mn concentration and an increase of the C/N ratio of the litter. Microarthropods could be of value for 
indicating a decline of soil quality with consequent loss of forest vitality. 

Key words: forest decline, soil acidification, manganese, nitrogen, C/N ratio, Oribatida, Platynothrus peltifer, 
Collembola, Sminthuridae. 
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